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^  SUMMARY 

Predictions  are  made  of  concentrations  of  NjO  produced  within  the  exhaust  of 

a  double-base  propellant  rocket  motor.  Typical  concentrations  produced  are 

. 

predicted  to  be  between ^10  and  10  molecule  ml  , 'jsome  five  orders  of  magnitude 
below  those  of  CO^  and  lower  than  the  characteristic  atmospheric  level  of  about  & 
1.3  x  10*^  molecule  ml  *. 
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In  INTRODUCTION 


Emissions  in  the  carbon  dioxide  Z*  (0,  0°,  I)  -  Z+  (0,  0°,  0)  band  and 

other,  overlapping,  CO^  bands  at  wavelengths  Ss  4.3  pm  are  a  prominent  infrared 

feature  of  rocket  exhaust  flames.  Nitrous  oxide  Z*  (0,  0,  1)  -  Z*  (0,  0,  0)  band 

emissions  also  occur  at  these  wavelengths,  and  absorption  of  both  N ^0  and  CO2 

radiations  results  from  the  presence  of  these  molecules  in  the  atmosphere.  Atmos- 
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pheric  concentrations  of  CO^  and  N^O  are  approximately  8  x  10  and  1.3  x  10 
molecule  ml  *  respectively*,  and  rocket  exhaust  COj  concentrations  may  be  calculated 
with  satisfactory  accuracy  by  the  methods  of  Ref.  2  and  3.  Exhaust  concentrations 
of  N2O  have  not  previously  been  calculated,  however.  This  memorandum  suomarizes 
results  of  a  calculation  of  NjO  production  in  a  particular  exhaust. 

2  METHOD  OF  CALCULATION 

Calculation  of  a  rocket  exhaust  structure  (the  distribution  in  space  of 
temperature,  species  concentrations,  etc.)  proceeds  in  three  stages.  First, 
conditions  in  the  chaaber  are  calculated4  on  the  basis  of  the  assumption  that 
chemical  equilibrium  is  achieved  at  the  high  temperatures  and  pressures  encountered 
therein.  Secondly,  the  no2zle  flow  is  treated  by  a  method  which  allows  for  non¬ 
equilibrium  chemical  effects  (a  simplified  one-dimensional  treatment4  was  used  in 
the  present  work),  and  the  nozzle  exit  conditions  are  obtained.  Thirdly,  the  exhaust 
structure  is  computed  from  knowledge  of  these  exit  conditions  and  such  external 
quantities  as  missile  velocity  and  altitude. 

\ 

The  third  stage  of  the  calculation  is  the  most  difficult.  The  technique  used 

2  3  5  *  ^ 

here  *  *  stems  from  coupling  together  a  two-equation  model  of  turbulence  and 

detailed  non-equilibrjtam  chemistry.  The  two  turbulence  variables  for  which  equations 

were  solved  in  this  work  were  the  turbulepce  kinetic  energy  k  (half  the  sum  of 

the  squares  of  RMS  fluctuating  parts  of  velocity  components)  and  W  ,  the  square 

of  a  frequency  characteristic  of  the  energy-containing  eddies. 

( 

The  governing  equations,  assumptions,  boundary  conditions  and  method  of  finite- 
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difference  solution  are  indicated  elsewhere  *  * 


With  the  notation  of  Ref.  2, 


the  turbulence  model  constants  used  here  took  the  following  values:  =  3.5, 

*  0.17,  Cj  -  1.48,  Cp  ■  0.09;  turbulent  Prandtl  and  Schmidt  numbers  Prfc  and 
Sc£  -  1.0,  corresponding  coefficients  for  k  and  W  -  0.86.  No  equation  was 
included  here  for  the  g  variable  of  Ref.  2.  The  set  of  reversible  chemical 
reactions  (with  rate  coefficients)  used  for  the  exhaust  calculation  is  shown  in 
Table  1.  Reactions  1-18  constitute  a  basic  exhaust  combustion  mechanism  (cf.  Ref.  5). 
Reactions  19-34  are  designed  to  account  for  production  of  NjO,  both  directly  from 
N^,  0^  and  0  and  indirectly  via  NO,  HN0,  N  and  NOj.  Rate  coefficients  for  reactions 
1-25  are  taken  from  Ref.  7;  those  for  reactions  26-34  derive  directly  from  Ref.  8. 
Thermochemical  data  for  all  species  stem  from  Ref.  9.  The  same  set  of  chemical 
data  was  used  for  the  nossle  calculation.  The  species  PbO,  SnO  and  Ca  were  treated, 
for  the  present  purpose,  as  chemically  inert,  PbO  and  SnO  being  assumed  to  be  present 
as  solid  particles  small  enough  to  follow  gas  streamlines. 

Exit  conditions  for  the  selected  (double-base  solid  propellant)  rocket  motor 

are  shown  in  Table  2.  The  calculation  was  performed  for  static  sea-level  conditions, 
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with  an  atmosphere  temperature  and  pressure  of  288  K  and  1.013  x  10  N  m  respect¬ 
ively  and  an  atmospheric  composition  of  Qn^  J  *  2.02  x  10^,  “  5.31  x  10^®, 

C  -1  *  ®*®®  *  1015  and  [other  species]  ■  0  molecule  ml  * .  (No  N2O  was  included 
in  the  atmospheric  composition  for  the  calculation  because  interest  was  centered  on 
production  of  this  species  in  the  exhaust  rather  than  its  entrainment  from  the 
atmosphere.)  The  exhaust  was  axisynane trie  and  contained  no  significant  region  of 
base  recirculation. 

3  RESULTS 

Selected  outputs  from  the  calculation  are  shown  in  Figs.  1-5.  Fig.  I  shows 

variations  of  [.NjO^l  with  radial  distance  r  from  the  exhaust  axis  for  specified 

axial  distances  x  downstream  of  the  nozzle  exit.  Peak  concentrations  of  NjO  are 
12  —1  ^ 

X  3  x  10  molecule  ml  ,  with  an  uncertainty  factor  probably  of  about  3.  Correspond¬ 
ing  radial  temperature  profiles  appear  in  Fig.  2.  Radial  profiles  of  COj  J  at 
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Reaction  mechanism  for  N20  calculations 


Reaction 


O  +  O  +  M  +  O^+M 
O  +  H+  M-*0H+M 

h+h  +  m-*h2+m 
H  ♦  OH  +  M  -*•  H20  +  H 
CO  +  0  +  M  -*■  C02  +  M 


OH  +  H„ 


H20  +  H 


0  +  H2  ->  OH  +  H 
H  +  02  ->  OH  +  0 
CO  +  OH  -»•  C02  +  H 
OH  +  OH  ->  H20  +  0 

co  +  o2  -*•  co2  +  0 


H  +  02  +  M  ■+  H02  +  M 
H  +  H02  ->-  OH  +  OH 
H  +  H02  -*•  H2  +  02 

H2  +  H02  H2°  +  0,1 
CO  +  ho2  ->  co2  +  OH 

o  ♦  ho2  -*■  oh  ♦  o2 

OH  +  H02  -*•  02  ♦  H20 


•  3 

Rate  Coefficient 

3  x  10~34  exp (900 /T) 

-29  -1 

1  x  10  *  T 

-30  -1 

3  x  10  T 

-25  -2 

1  x  io  "  T  4 

7  x  IO-33  exp(-2200/T) 


1.9  x  IO-15  T1*3  exp(-1825/T) 
3  x  io"14  T  exp(-4480/T) 

2.4  x  io-10  exp (-82 50 /T) 

2.8  x  io"17  T1*3  exp(330/T) 

1  x  io"11  exp(-550/T) 

4.2  x  io  ^2  exp(-24000/T) 


2  x  io"32  exp (500 /T) 

4  x  io"10  exp(-950/T) 

4  x  io'11  exp(-350/T) 

1  x  io“12  exp(-9400/T) 

2.5  x  10  ^  exp(-11900/T) 
8  x  IO-11  exp (-500 /T) 

5  x  io-11 


NO  +  NO 
N20  ♦  H 


n2o  ♦  0 

N2  +  OH 


NO  ♦  H  ♦  M  HNO  ♦  M 
HNO  ♦  H  -*■  NO  ♦  H~ 


HNO  ♦  OH 

n2  ♦  0 


NO  ♦  HjO 
NO  ♦  N 


N  ♦  0,  -*■  NO  +  0 


N2  ♦  °2  ♦  N20  ♦  0 
N2  ♦  0  ♦  M  ■*>  N20  ♦  M 

no2  ♦  o  -*■  no  ♦  o2 

N02  ♦  H  -*>  NO  ♦  OH 
n  +  n  +  m*n2  +  m 

N  ♦0*K-»IW  +  M 


NO  ♦  0  ♦  H 


«♦  N02  ♦  M 


N  ♦  ON  •*>  NO  ♦  H 


HNO  ♦  WO 


♦  m2o  ♦  h2o 


2.2  x  io‘12  exp(-32100/T) 

1.3  x  IO*10  exp (-7600 /T) 

5  x  10“32  axp(300/T) 

8  x  IO* 12 

6  x  lo"U 

1.3  x  io“10  exp (-38000 /T) 
1.1  x  io**14  T  exp(-3150/T) 
1.1  x  IO"10  exp(-55300/T) 

4  x  10~35  exp (- 10660 /T) 

1.7  x  IO-11  exp(-300/T) 

5.8  x  10“10  exp(-740/T) 

8.3  x  10~34  exp(500/T) 

.aw  ,  -31  _-0. 5 

1.8  x  10  T 

3  x  io”33  exp(940/T) 

7  x  10"U 

1  x  io“12  axp(-2000/T) 


TABLE  2 

Rocket  Motor  Exit  Plane  Conditions 


1  Quantity  and  Units 

Value 

Species  concentration ,  molecule  ml 

c°2 

2.2  x  1018 

N 

1.8  x  1018 

H,0 

3.3  x  1018 

CO 

4.7  x  1018 

H2 

1.4  x  io18 

H 

8.7  x  IO15 

* 

OH 

4.7  x  IO14 

°2 

12 

7.0  x  IO14 

0 

12 

2.4  x  10 

H02 

9.3  x  IO9 

n2o 

9.4  x  IO9 

N 

4.3  x  IO9 

NO 

9.1  x  IO14 

NO, 

8.0  x  IO7 

4 

HNO 

13 

2.6  x  10 

' 

PbO 

3.3  x  IO16 

SnO 

1.8  x  1016 

Ca 

5.0  x  IO15 

i 

1  Jet  velocity,  km  8 

2.2 

■  -2 

5 

■  Pressure,  N  m 

2.7  x  10 

1  Temperature,  K 

1470 

1  Nozzle  radius,  cm 

5.74 

■  2-2 

1  Turbulence  kinetic  energy  k  ,  m  s 

7  x  IO4 

W  ,  • 

Mass  flow  rate,  kg  s' 


2.1  x  ior 
13.3 


wmmmmmmmmmmmmmim' 


the  same  axial  stations  are  given  in  Fig.  3;  comparison  between  Fig.  1  and  Fig.  3 
reveals  that  typically  exceeds  L^oll  by  about  five  orders  of  magni¬ 

tude  in  the  hotter  parts  of  the  exhaust.  Axial  profiles  of  N20^|  an<*  tempera¬ 
ture  are  shown  in  Figs.  4  and  5  respectively:  both  the  shock  structure  close 
to  the  nozzle  exit  and  the  boost  in  temperature  resulting  from  secondary  combustion 
of  CO  and  H2  are  clearly  apparent. 

Two  subsidiary  computations  were  performed.  In  the  first,  the  sensitivity  of 
results  to  Prfc  and  Sct  was  illustrated  via  setting  of  these  numbers  to  0.7 
rather  than  1.0,  all  other  inputs  remaining  unchanged.  The  distribution  of  tempera¬ 
ture  (Fig.  5)  and  of  £n20"2  (Fig.  6;  compare  Fig.  1)  were  not  greatly  affected, 
although  the  exhaust  plume  became  slightly  shorter  and  narrower.  In  the  second, 
values  of  0.7  for  Prt  and  Scfc  were  retained  but  all  chemical  reactions  of 
nitrogen-containing  species  except  26  and  27  were  omitted.  This  change,  designed 
to  isolate  effects  on  °f  two  reactions  regarded  at  the  outset  as  likely 

to  play  dominant  direct  roles  in  determining  this  concentration,  caused  the  predicted 
maximum  concentration  to  rise  by  a  factor  of  2-3  (Fig.  7).  Reactions  other  than 
26  and  27  evidently  do  play  significant  parts  in  determining  \~  ^0  J  . 

It  is  worth  emphasising  that  N20  concentrations  produced  in  the  exhaust  are 
well  away  from  those  consistent  with  local  thermochemical  equilibrium.  The  last 
three  columns  of  Table  3  show:  (a)  concentrations  F  27  calculated  on  the 
assumption  that  reaction  27  is  locally  balanced;  (b)  concentrations  Q^oll  2$ 
calculated  on  the  aasumption  that  reaction  26  is  locally  balanced;  and  (c)  concentra¬ 
tions  f N20~|  derived  from  the  (standard  inputs)  exhaust  structure  calculation. 
l~N2°II  differs  markedly  from  both  CN2°U27  and  £n2°“| 26  because  non- 
equilibriiaa  effects  (especially  those  affecting  [oj  )  are  critically  important. 


v 
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4  DISCUSSION  AND  CONCLUSIONS 


12 

The  maximum  concentration  of  ^0  produced  in  this  exhaust  is  %  3  x  10 

-1  13  -1 

molecule  ml  .  This  is  lower  than  the  level  of  1.3  x  10  molecule  ml 

characteristically  present  on  the  atmosphere.  The  exhaust  NjO  concentrations, 
even  if  raised  somewhat  via  allowance  for  entrainment  of  atmospheric  NjO  ,  are 
some  five  orders  of  magnitude  below  those  of  CO2  .  It  is  likely  that  the  main 
effect  of  N2O  on  IR  radiations  at  wavelengths  *  4.3  pm  will  be  to  cause  atmos¬ 
pheric  attenuation  effects  over  long  path  lengths  at  certain  specific  wavelengths 
rather  than  to  produce  significant  contributions  to  (exhaust)  source  radiation. 
Such  a  conclusion  may  reasonably  be  expected  to  hold  for  a  wide  range  of  rocket 
motors . 
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Predictions  are  made  of  concentrations  of  NjO  produced  within  the  exhaust  of  a 
double-base  propellant  rocket  motor.  Tynical  concentrations  produced  are  predicted 
to  be  between  10u  and  10 13  molecule  ml-1,  some  five  orders  of  magnitude  below  those 
of  CO,  and  lower  than  the  characteristic  atmoapharic  laval  of  about  1.3  x  lO1^  I 
molecule  ml"1.  1 


